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Abstract. The two major discrete time formulations for quantum walks, coined 

and scattering, are unitarily equivalent for arbitrary position dependent transition 

C^ . amplitudes and any topology (PRA 80, 052301 (2009)). Although the proof explicit 

^ , ' describes the mapping obtention, its high technicality may hinder relevant physical 

aspects involved in the equivalence. Discussing concrete examples - the most general 
constructions for the line, square and honeycomb lattices - here we unveil the 
similarities and differences of these two versions of quantum walks. We moreover show 

f~^ . how to derive the dynamics of one from the other by means of proper projections. 

^^ I We perform calculations for different probability amplitudes like, Hadamard, Grover, 

Discrete Fourier Transform and the uncommon in the area (but interesting) Discrete 
Hartley Transform, comparing the evolutions. Our study illustrates the models 
interplay, an important issue for implementations and applications of such systems. 
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1. Introduction 

Since their introduction [HE], [3], quantum walks (QWs), a paradigmatic and relatively 
simple class of quantum systems |3] , have found many applications in different areas of 
science (see, e.g., the recent review in [5j). At a first glance, such usefulness may be 
attributed to QWs distinct possible formulations. For instance, although in all cases 
the dynamics take place in discrete spatial structures, graphs (or lattices), time can be 
either a continuous [6j or a discrete variable. Moreover, in discrete time versions, the 
evolution can be dictated by inner states called "coins" |7] (coined QWs, CQWs) or 
resulting from scattering-like processes [HI El ED] (scattering QWs, SQWs). 

Regarding usages for QWs, some aspects of the general searching problem [11] 
seem to be very appropriate for SQWs [12], including identification and comparison 
of parts of a graph [13], [Tl]. Also, SQWs are frequently considered in the study of 
scattering in semi- infinite graphs [9] [TOl |15], a relevant configuration to implement 
universal quantum computation through QWs [I6]. On the other hand, the literature 
proposes a much larger number of applications for CQWs [5] (TT] [HI [19]. Examples 
are: transport in biological systems [20] (in a continuous time context, see [21]); Bose- 
Einstein condensates redistribution [22]; quantum phase transition in optical lattices 
[23j : decoherence processes [21]; and even quantum games [25]. In particular, CQWs 
are widely discussed in quantum computation as a tool for the development of quantum 
algorithms [26]. 

As it concerns implementations, distinct experimental setups based on distinct 
physical phenomena can be used to actually build QWs. For CQWs, protocols 
based on trapped atoms in optical lattices [27J, quantum dots [28], photons orbital 
angular momentum [29], and QED cavities [30], to cite just a few, have been devised. 
Concrete lab realizations were constructed with waveguide lattices [31] (refers to [32] for 
waveguides in the continuous time case), passive optical elements [33], and liquid-state 
nuclear magnetic resonance [31], again only mentioning few examples (an overview is 
given in [35]). By their turn, SQWs could be associated to optical networks [36] [37] 
and eventually may be ensemble just with linear optical elements [121 [13], similarly as 
done in [38] to fabricate a quantum version of Galton's quincunx (a classical mechanical 
machine which at certain locations "chooses" - with 50%:50% probability - between 
two directions to go for a traveling ball). 

The above mentioned models distinctions for applications (and in a less extend for 
implementations) are, however, due to practical and operation instead to fundamental 
reasons. Indeed, continuous time and coined QWs are closely related, since they lead 
to a same dynamics in proper limits [391 SO] (although the limits may not be so direct 
to achieve [IQ]). Furthermore, it has been rigorously proven that CQWs and SQWs 
are unconditionally (i.e., for arbitrary topologies and spatially dependent transition 
probabilities) unitary equivalent [il] . 

These kindred relations open important perspectives in the employment of QWs, 
specially for the discrete time versions. They should be equally appropriate in any 
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application and in principle (see comments in Sec. 4), a same physical implementation 
would be capable to simulate both CQWs and SQWs. But to benefit from such 
connections, it is necessary a clear understanding of similarities and differences of the two 
system versions as well as the correct way to obtain one from another. The proof in |41j . 
having as the main goal to formally demonstrate their equivalence, relies on general and 
very technical arguments, making difficult to fully appreciate the conceptual features 
associated to the QWs reciprocities. Thus, knowing that CQWs and SQWs are always 
equivalent, we can focus on representative cases, discussing in details the interplay 
between the discrete time formulations. Hence, here we illustrate that different QWs 
are in fact closely tied and a specific model choice might be much more a matter of 
practicality than of unfeasibility of the corresponding sibling construction. 

The work is organized as the following. From a parallel with usual classical random 
walks, in Sec. 2 we discuss for the ID case the main characteristics of CQWs and 
SQWs, also addressing their mapping. We do so assuming the most general situation 
of complete arbitrary position dependent transition probabilities at the ID lattice sites. 
In Sec. 3, the same type of analysis is carried out for two 2D lattices, the square and 
the honeycomb. The models are formulated in details and simulations, exemplifying 
their distinctions for the time evolution along the graphs structures, are presented. 
Importantly, we show how the spatial probabilities distributions for one model can be 
obtained from the other through proper projections. Finally, results for different coin 
(and the equivalent scattering) matrices operators are compared. It includes the more 
common Hadamard, Grove and Discrete Fourier Transform, as well as the not so usual 
(but interesting) Discrete Hartley Transform and few others. A short conclusion is 
drawn in Sec. 4. 

2. Discrete quantum walks on the line: two different views 

To better understand the main ideas underlying the construction of the existing discrete 
QWs models, we first recall an interesting way to view the classical case on the line (ID). 

Consider a classical walker, starting in x = at t = 0, that takes steps of fixed 
length L and moves with constant speed v. Its simple dynamics is described as the 
following. Each time the walker reaches a position x{t = nr) = x„ = ±jL (j = 0, 1, . . .) 
for n = 0, 1, . . . and r = L/v, it randomly (and instantaneously) chooses a new direction 
to go, either to the right (a = +1), with probability p, or to the left {a = —1), with 
probability 1 — p. So, this ID continuous system can be characterized by two processes 
taking place in an "effective lattice": (i) one purely stochastic (choosing directions), 
occurring at the "lattice sites" ±jL; and (ii) other purely deterministic (ballistic motion, 
with x(t) = Xn-i + av{t — tn-i) for (n — l)r < t < nr), occurring along the "lattice 
bonds". Such view is schematically represented in Fig. 1. 

The construction of a discrete time QW closely follows the previous picture. 
Essentially, there are two possible implementations, based on what we assume as the 
primary process, either (i) or (ii) above, to describe the quantum states. Indeed, in the 
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Figure 1. (a) A classical random walk faced as alternate stochastic and deterministic 
processes along a ID lattice, (b) Example of a possible path for t = 6r, whose 
probability is p^(l — p)^. Here, p [1 — p] is the probability to go to right [left]. 
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Figure 2. For QWs in ID, the associated "Hilbert lattice" , but which is not necessarily 
a spatial structure once the states do not need to be position eigenvectors. The states 
are defined (a) on the sites, coined, and (b) along the bonds, scattering, formulations. 
For SQWs, details in (b) show the scattering quantum amplitudes defined on each site. 



first case (associated to CQWs) the quantum states describe the system at the lattice 
site positions: classically, the locations where it is made a probabilistic choice about 
the next step direction. In the second case (associated to SQWs) the quantum states 
are defined on the bonds: classically, corresponding to the deterministic locomotion 
along the lattice. This is an important distinction since each model has a different state 
representation with a different possible interpretation. 

One does not need to relate any "lattice" structure to QWs. Nevertheless, it is 
instructive to go on with such analogy and to view quantum walks in ID as a dynamics 
defined on a "Hilbert lattice" [12], depicted in Fig. 2. In this way, we can associate 
the classical step length L to the characteristic unit distance Aj = 1 between two 
consecutive "sites" of the Hilbert lattice. Moreover, in both cases the discrete time 
evolution is given by a single step unitary operator ?7, acting on states belonging to this 
Hilbert (lattice) space, so that |^(n-h 1)) = U\^{n)). 
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2.1. The coined quantum walk model 

Quantum states are defined on the sites j, Fig. 2 (a). Then, we consider the basis states 
{|j)}, which span the Hilbert subspace Tip. Moreover, to any j we associate an inner 
two- level Hilbert subspace Tic, whose "coin" states are |±). They give the correct unitary 
"stochastic" character for the problem [21 03] . Also, their role is somehow similar to that 
of choosing directions in the classical case: once \j) represents only spatial "location" 
(sites), then "direction" must be played by the auxiliary coin states. 

The appropriate orthonormal basis to describe the entire space Ti = Tip ® Tic 
(L2(Z) ® L2(Z2)) is {\j) ® \a)}, with a = ± and j = 0, ±1, ±2, . . .. Defining a shift 
operator as S\j) = \j + 1) (and S^j) = |j — 1)), we write 

U,= {S0 |+)(+| + 5t ® |-)(-|) (t \j){j\ ® C^A ■ (1) 

Uc is unitary if the coin-flip operators C*^-^-* 's are also unitary. There are many possible 
choices for the C's [44], a common one being the Hadamard's, for which V j, C^^' = H2 
with H2\a) = {\ — a) — a\a))/\/2. However, the most general case is to consider the 
matrix elements ((t"|C*^-'^|(t') = c)^,l^, and to assure unitarity just by imposing 

I Ji) |2 , I (i) |2 _ I (i) |2 , I (j) |2 _ 1 I (i) |2 _ I (i) |2 

I t^_^ _j_ I ~r I ^ |- I — I *-■— — I ~r I L--|- _ I — -L 5 I L-^ _ I — I L._ _j_ I , 

cJ'V[c^'V]* + c?-[c^l]* = 0. (2) 

Thus, we have 

Uc\j) ® W) = c^£ \j + a)0 \a) + c^ll„ \j-a)®\-a). (3) 

One of the most important applications of CQWs is in quantum computation. So, 
it is usual to define the above matrix elements in a way easy to associate to quantum 
gate operators. It explains, e.g., the frequent choice in the literature of the Hadamard's 
coin (fundamental to manipulate q-bits). This practice, however, leads to constructions 
of coin matrices in a format not common among physicists, but popular in mathematics 
and computer science (as nicely explained and exemplified in Ref. |15], e.g., p. xxiii). 
Indeed, the states |±) are written as the column matrices 



)= . l+)= . (4) 





opposite to the usual spin-up and spin-down convention. Following this representation, 
the most general form for the coin operator, obeying to Eq. ([2]) and having the 
Hadamard as a particular case is (see also [IHl HTJ HE]) 

r(i)=p^^r,^l (^ 6^P[^^i]cos[^,-] exp[zg] sin[^,-] A 

^^^^^ 1^ expHO]sin[^,] -expHej]cos[0,] y • ^' 

with < IjXj^ij^^j < 27r. Then, Eq. (|3]) reads 

Uc\j) ® |cr) = - (jexp[i(7j - a^j)] cos[^j] \i + o) ®\+ a) 

+ exp[+z(7j + aCj)] sin[^j] \j - cr) ®\- a). (6) 
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2.2. The scattering quantum walk model 

For our purposes, it is easier to formulate the scattering model in a slightly different 
(but completely equivalent) way than usually done in the literature [10]. Along each 
bond, joining two consecutive sites, say j — 1 and j (see Fig. 2 (b)), we assume two 
possible states | + l,j) and | — l,j — 1). Hence, contrasting to the coined model, now 
the quantum number a = ±1 is associated to a "direction" along a bond in the Hilbert 
lattice. The full Hilbert space "H is no longer a direct product of two subspaces. Each 
basis element |cr, j) belongs to L^(Z x Z2) and satisfies to {j',o''\a,j) = 6jij6„io-- 
Defining the operators T and R by 

T\a,j) = t(/) \a,j + a), R\a,j) = r^^^ \-a,j- a), (7) 

and 

Tt|a, j) = t^'^)* \a,j - a), R^a,j) = riV^^^* | -a,j- a), (8) 

the one step time evolution is simply 

Us = T + R, (9) 

so that 

Us\a,j) = t^^ lor, J + or) + r^^ I - a, J - a). (10) 

The unitarity of Us implies [19] 

|t0-)|2 + 1^0-) |2 = 1^ 1^0-) |2 = |/.-)p, rStt^')* + r(^')*t(^-) = 0, (11) 

which are exactly the relations satisfied by the reflection and transmission amplitudes 
in a quantum scattering problem in ID ^91 [50] , resulting from the unitarity of the S 
scattering matrix. 

Equation (ITTl) automatically holds if for any j (with < pj < 1 and < Aj, (pj, ipj < 
27r) 

t^^) = exp[zAj] (a/1-Pj exp[za-0j]) , 



r 



(i) 



exp[iAj] {cr ^JJTj exp[iav9j]) . (12) 

If \j = A for all j, without loss of generality we can set A = 0. We should mention 
that the expressions in Eq. flT2|) are not the only possibility [12]. There is a certain 
arbitrariness in signals convention. The present choice, however, has a direct physical 
motivation. For time-reverse invariant systems t^ = t_ ^0]. Then, if we also assume 
real amplitudes, the signals for the r's must be opposite, exactly the case in Eq. ( !T2l) . 
Furthermore, for j-independent scattering coefficients, the above (with (pj = cpj = 0) 
are the relations adopted in the original work introducing SQW models [8]. 

2.3. Obtaining the probabilities 

At this point we should emphasize that the above constructions are more general than 
simply to quantize the dynamics of a common classical random walk. Here, by common 
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we mean that each time a new direction needs to be chosen, we use the same probabihties 
{p and 1 — p) to decide between right and left. Quantum mechanically, to allow the 
parameters to depend on j (cf., Eqs. ([5]) and ( TT2l) ) implies that we explicit assume 
position dependent probability amplitudes. Obviously, by making such parameters j- 
independent, we recover the relation with the usual classical case. 

Now, suppose we shall determine which is the probability P^^\n) to be in the 
"position" state j (which means a site (bond) state in the coin (scattering) model) at 
time n, regardless the value of the coin (direction) quantum number a. We define then 

^?'^ = ij)(ji®Ei^)(^i' n^^'^ = Ei^'-^' + ^)(^' + ^'^i' (13) 

(T (J 

respectively, the coin and scattering position projector operators. The sought probability 
is the expected value 

P^\n) = {^{n)\V^\^{n)), (14) 

for V'^^^ one of the expressions in Eq. (fT3ll . 

Projection is thus an essential ingredient in defining a QW model, which differs 
from classical walk systems by the typical interference effects due to the Eq. flT4|) (see, 
e.g., Ref. [51]). 

2.J^. The unitary equivalence of the two ID quantum walk models 

Finally, in the present ID topology the unitary equivalence between CQWs and SQWs 
is straightforwardly established. 

First, note the one-to-one correspondence between their full Hilbert spaces, L^(Z)® 
L^(Z2) = L^(Z X Z2). So, it follows directly the existence of an isomorphic unitary 
operator [H] -E : "H — )■ "H, given by 

E W,j) = \j) ® \a), E^ \3) ® \a) = |a, j). (15) 

Second, (i) inspecting Eqs. ([3]) and (fTOj) : (ii) considering the most general conditions 
for the models quantum coefficients, represented by Eqs. ([2]) and (iTTl) (observe that Eqs. 
([5]) and flT2l) . more usual in the literature, already obey particular conventions); and 
(iii) taking into account Eq. ( TT5l) : one finds that the following coefficients play complete 
similar roles 

c'^l^t^^^ and c^H.^r^jl (16) 

Hence, setting Cal = tj and c_l„ = rj for all j, we have in both models exactly the 
same probability amplitudes for their time evolutions. 

Third, from Eqs. (I3l), ( TTOj) . ( ITSl) and assuming the above equalities, the resulting 
dynamics are unitary equivalent once 

Us = E^U,E. (17) 

As it is known [8j, even when different QW models are unitary equivalent, the 
associated probabilities - obtained through direct projections - can be distinct. This 
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is so because we are choosing contrasting physical representations, sites (coined) and 
bonds (scattering), to describe the problem. The specific states which characterize one 
of these "spatial" configurations are not akin to the states for the other. For instance, 
the two states corresponding to the site j in CQWs, |j) ® |cr = ±), are mapped to states 
at different bonds in SQWs. Therefore, the probability to be in a unique site is not 
equal to the probability to be in a unique bond (cf., Eqs. flT3|) - flT^ ). Mathematically, it 
is related to the fact that considering Eqs. flT3|) and flTSj) . we get Vs 7^ E'^Vc E and 
Vc 7^ EVs E^ (inequalities holding true in any topology given the results in [H]). 
Thus, for |\E')^ and \^)c = E 1"^)^, we have from Eq. (^ that generally P]^^ 7^ p!^^K 

Nevertheless, due to the QWs unitary equivalence we can recover the probabilities 
from each other version by constructing proper cross projector operators. Indeed, 
defining 

pO') 1^ = E^ V[^^ E, Vi^^\^ = E V\^^ E\ (18) 

one readily obtains the probabilities for the coined (scattering) model by applying Ps [^ 
(Pc Is) to the state \^{n))s (|\l'(n))c), evolved according to the scattering (coined) 
formulation. 

The numerical examples for 2D lattices in Section 3 will clearly illustrate all the 
above observations. 

3. More general topologies: two illustrative examples 



As already mentioned, it has been proven in Ref. [H] that CQWs and SQWs are always 
related by an unitary transformation. In the previous Section, profiting from its relative 
simplicity (although considering the most general situation, hence extending the results 
of Ref. [8]), we have explicit illustrated so for the ID case. 

However, the benefits in being able to map CQWs and SQWs become really evident 
when one considers more complex topologies. Thus, next we present a detailed analysis 
of the correspondence between CQWs and SQWs in two particular, but very instructive, 
2D examples, square and honeycomb lattices. 

3.1. Quantum walks on a square lattice 

The square lattice, represented in Fig. 3, is the natural 2D extension of the ID topology. 
In the following we describe in such a case the coined and scattering QW formulations, 
as well as their unitary equivalence. 

3.1.1. The coined formulation 

In Fig. 3 (a) we schematically represent the coined Hilbert lattice, whose basis 
states are {\j) ® \k) ® |cr)}, with a = 1,2,3,4 and j,k = 0,±1,±2,.... For notation 
simplicity, we do not write the "spatial" states in the form \j)x^ \k)y (and operators as 
Ax ® By) to distinguish between the two distinct ID dimensions. Also, we assume the 
natural ordering convention {A® B ^ C) \a) ® |/3) (g) I7) = A\a) ® B\l3) (g) CI7). 
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Figure 3. A QW on a square lattice, (a) At each site (j, fc), the arrows schematically 
represent the four inner states, a = 1,2,3,4, of the coined formulation, (b) For the 
scattering formulation, it is shown the four "incoming" states propagating to the lattice 
sites marked in (a) . The two possible states defined on the bond common to the regions 
A and B (A and C) are |++, (j, k)) and | , (j-l, k)) (|-+, (j, k)) and |+-, (j, fc-l)>). 
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The unitarity of C, directly extending the relations in Eq. ([2]), leads to 
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Figure 4. Schematics of tlie one step time evolution for tlie coined model in a square 
lattice in the natural (a) and the diagonal (b) cases. The system, initially at the state 
represented by the site x , in a first step can reach the states represented by the sites 
■ . In the second step, the states which can be visited are those represented by the 
sites • and x . The diagonal version (dotted lattices) can be thought as the natural 
one rotated by 7r/4 and rescaled by a factor v^- 



As the one step time evolution, there are two common possibihties. The first is the 
so called natural choice, since single steps displacements follow the natural topology of 
the lattice [52], i.e., from a given site representing a state, it can go either to right, left, 
up or down [H |S3] (Fig. 4 (a)). The second is the diagonal version [SI]- In this case, 
single steps are those in which "moves" take place along the diagonals of the Hilbert 
lattice (Fig. 4 (b)). 

Let us define (for 5* as in Sec. 2-1) 

AT = 5 ® I ® 1 1) (1| + 5^ ® I ® |3) (3| + I ® 5 ® |2) (2| + 1 ® 5^ ® |4) (4|, (22) 
with I the identity operator. The time evolution for the natural version, an extension 
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of Eq. (HD, yields 

Uc = Nx (5^|j)(j|®|fc)(fc|®C(^''=)). (23) 

For the diagonal version, we first define the operator 

D = I®5®|1)(1| + I®5^®|3)(3| + 5^®I®|2)(2| + ^®I®|4)(4|.(24) 
So, the time evolution is given by 

'--'c— diagonal ^ ^-^ ^ c- V'^"-'/ 

Equation (1251) . with D as in Eq. ( 12^ . has a very simple interpretation (if we maintain 
that Eq. fl25|) represents a single step). Indeed, the diagonal case is essentially the 
natural one rotated by 7r/4 (anti-clockwise) in the Hilbert lattice. Moreover, quantities 
associated to any norm (length scale) defined on the lattice, e.g., diffusiveness, should 
be re-scaled by a factor a/2, which is just the ratio between the diagonal and side of 
the unitary cell. Fig. 4. In all our further analysis, we will consider only the natural 
version, observing that the diagonal would easily follow from Eq. (1251) . 
Finally, for the natural coined model the projector operator is 

4 

^P^ = |j)(j|®|A:)(A:|®5^k)(a|. (26) 

o-=l 

3.1.2. The scattering model 

For (Tj;, ay = ±1 (± for short) and j, k = 0, ±1, ±2, ..., the scattering model basis 
states are {\crxcry, (j, k))}, Fig. 3 (b). The time evolution is again Us = T + R, with 

(27) 

The corresponding actions of T"^ and R^ are a straightforward generalization of Eq. (|8]), 
in a format similar to Eq. ( 1271) . Defining 

/ ^(i.fc) /J,h) ^{j,k) ^(j,k) \ 



T axay,{j,k)) = 




Raxay,{j,k)) = 


— 0"x —Cy, Gx O'y 



Jj,k) .(J,k) (j,k) (j,k) 

I J L __ t-J L J L t'J U J t-J U _ 



Ahk) Aj,k) Aj,k) {3,k) 

''+-, ''+-,++ ''+-,+- '+-,-+ 

V Ahk) .U,k) ij,k) Aj,k) , 

\ ''-+, ''-+,++ '-+,+- ''-+,-+ / 



(2J 



the unitarity of Us implies in the unitarity of the scattering matrices F*^-'''^'' (once more a 
direct extension of previous results, in this case of Eq. (TTTj) ). So, the elements of F*^-'''^'' 
must satisfy to relations completely analog to those in Eq. ( |2T1) . 
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For the projector operators, we need to distinguish between the horizontal (a^ x ay 



+) and the vertical (cr^ x a^ 



bonds. Hence, we define 



'pU,^;^)=J2W^^^,{j + 



cr^ 



k)){ij + 



cr^ 



k),a^a^ 



vi^''-'-'> =Y.\-^y ^^' (3' ^ + ^))(0' k 



S-i- 



, ayay\. 



(29) 



3.1.3. The mapping 

To map the two models, we first make the following identification between the 
scattering directions and the inner coin quantum numbers {a^CTy ^ a): 

+ + ^ 1, -+ o 2, o 3, +- ^ 4, (30) 

which can be cast as a = (5 — (2 + ax)ciy)/2. 

Using the procedure in the last Section (or likewise, the rigorous construction in 
Ref. 5^), i.e., to associate the scattering states "incoming" to a certain site (see Fig. 
3 (b)) with the coin states at that site, we can set E as (taking into account Eq. fl5Ul) ) 



5 (2 
E\a,ay{j,k)) = \j)^\k)®\--^ 



cTxjcry 



)• 



(31) 



Furthermore, assuming the coefficients in Eq. (I20l) equal to those in Eq. (I28l) (which is 
consistent with the relation in Eq. fl30|) ). again we find that the two models are unitary 
equivalent, since for Eqs. ( !23l) . ( l27l) and ( l3Tl) . the relation in Eq. f ITTj) holds. 



Finally, to obtain the cross projector operators, we define Vs \ and Vc L as in 



U,k) I 



>(i,fc) I 



>(i,fc) 



ij,k) 



Eq. ([HD, for E given by Eq. §2), V^' ' by Eq. ([26]), and Vr ' by the appropriate 
pF''^)'^ in Eq. 



3.1.4. Examples 

To illustrate the above general constructions, we analyze some particular cases for 
the probability amphtudes (see, for instance, Refs. [TJ [521 [53]). For simplicity, we 
suppose all the coin (and therefore the scattering) matrices to be independent on the 
quantum numbers j and k. We choose the following four coin operators (whose corre- 
sponding scattering matrices are written in exactly the same form): 



(a) x-y decoupled Hadamard 



C 



H21SH2 




1 

7! 



/ +1 +1 \ 
+1-10 
0+1+1 

y +1 -1 / 



(b) Full 4x4 Hadamard 

/ +1 +1 +1 +1 \ 



C 



Hi 



-1 -1 +1 -1 

4 +1 -1 -1 
-1 -1 -1 +1 



(32) 



(33) 
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fc) 4 X 4 Grover 



C, 



G4 



/-I 

+1 
+1 

V+1 



+1 
-1 
+1 
+1 



+1 
+1 
-1 
+1 



+l\ 

+1 
+1 

-1/ 



(34) 



(35) 



(d) 4 X 4 Discrete Fourier Transform (DFT) 

/ +1 +1 +1 +1 \ 

^ _ 1 +1 +i -1 -i 
Cnrn-^ +1 _i +i _i 

y +1 -i -1 +i ) 

We compute t/" |\l/(0)), with n = 20, for both QW models and the above matrices. 
For each (j, A;) we use proper projectors to calculate the probability P'^^'^^{n = 20) = 
(^(20)|p(^''=)|^(20)). As 1^(0)), we take 

l*(0)). = ^[| + + (0,0)) + z|--(0,0)) + |- + (0,0)) + z| + -(0,0))], 

1 



|*(0))e = ^|^(0)), 



|0)®|0)® |l) + z|3) + |2) + z|4) 



(36) 



The p(J''^)'s are displayed as 3D density plots. So, in each graph the heights in the 
z-axis correspond to the probabilities values. Moreover, here and for the honeycomb 
lattice in Sec. 3.2, each pair of labels j and k in the x-y plane indicates the [j, k) lattice 
site spatial location. Thus, for the coined formulation, the p(-?''^)'s are marked just over 
the sites. For the scattering formulation, since the states are defined along the bonds, 
we mark the P's exactly over the middle points of the corresponding bonds. Hence, we 
show the probability patterns of each QW version in its own state representation, but 
in a way which makes easy to qualitatively compare the two models. 

The results for the decoupled Hadamard, full Hadamard, Grove, and DFT operators 
are presented, respectively, in Figs. (5) to (8). The plots are organized as the following. 
For the coined formulation. Figs. 5-8 (a) [(b)] show the probabilities obtained from the 



projectors in Eq. fl2^ [the projectors V, 



ihk) I 



On the other hand, for the scattering 



formulation. Figs. 5-8 (c) [(d)] show the p(-?''^)'s from the projectors in Eq. (129!) [the 
projectors Vs' |J. 

By inspecting Figs. (5)- (8), some points become evident, (i) Due to decoupled 
form of the transition amplitude matrix, Eq. fl32|) . in Fig. 5 we actually have two 
independent one dimensional evolutions. Hence, each direction is a good example of the 
ID results in Sec. 2. (ii) Despite the fact the two QW versions are completely equivalent 
by unitary transformations, they lead to distinct probability patterns (compare (a) or 
(d) with (b) or (c)). Indeed, as already emphasized (for instance, the discussion just 
after Eq. flT7|) in Sec. 2.4), the one-to-one correspondence between the two formulations 
pinpoint a dynamical similarity. However, the states in each case represent different 
(although correlated) aspects of the Hilbert lattice, thus the distributions of the P's 
do not need to coincide, (iii) Moreover, such patterns display richer structures in the 
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Figure 5. Quantum walks in tlie square lattice considering the decoupled Hadamard 
transition coefficients, Eq. (|32|) . The density plots represent the probabilities to be 
in the difi^erent states (as defined by the projector operators) after n = 20 time steps 
for 1^(0)) given in the main text. The results are for the: (a) coined, (b) scattering 
obtained from the coined, (c) scattering, and (d) coined obtained form the scattering, 
models. The graphs (a) and (d) and (b) and (c) are complete identical. 



scattering case. To understand so, assume the simpler ID lattice, for which the initial 
state |+,0) (£'|+,0) = |0) ® |+)) evolves one time step. For the SQW (CQW), we 
get 1^(1)) = rf I-, -1) + tf |+, +1) (1^(1)) = c% I - 1) ® I-) + cf+ I + 1) ® 1+)). 
Considering the Hilbert lattice picture, note that for the former we have two "neighbor" 
bond states in the expansion of |\E'(1)), those "attached" to the site 0. On the other 
hand, for the latter the site states composing |\E'(1)) are j = —1 and j = +1, but not 
j = 0. It illustrates a very typical situation in any topology, namely, the SQW dynamics 
tends to excite "contiguous" spatial (bond) basis states, whereas CQWs may skip some 
"successive" (site) basis states. Thus, interference |5l] is usually more recurrent in 
SQWs than in CQWs (due to this difference in the spreading of |^)), explaining the 
behavior observed in the plots, (iv) Nevertheless, we can recover the p(J''^)'s from each 
other model by means of the cross projectors (e.g., SQWs in (b) from CQWs in (a) 
and CQWs in (d) from SQWs in (c)), since the information about one model is always 
encoded in the other. 

We finally mention we have analyzed the diagonal CQW and its corresponding 
scattering version. We have obtained one formulation P's from the other by correctly 
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P 0.5l 



Figure 6. The same as in Fig. 5, but for tlie 4x4 Hadamard, Eq. ([331 




Figure 7. The same as in Fig. 5, but for the 4x4 Grover, Eq. ([M|) . 
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Figure 8. The same as in Fig. 5, but for tlie 4 x 4 DFT, Eq. (|35l) . 

defining the cross projectors. Furtliermore, as expected tlie probabilities plots (not 
shown) are exactly the ones here, just rotated by 7r/4 and rescaled by a factor a/2. 



3.2. Quantum walks on a honeycomb lattice 

Our last example is a QW defined on a honeycomb lattice, whose structure is depicted 
in Fig. 9 (a). It is far more involving than the previous square case mainly because 
now the coordination number, equal to 3, is odd |17] (for coordination 3 in ID see, 
e.g., Ref. [55] )• We note that although this topology was recently investigated in the 
context of continuous time models [56], and few discrete time implementations do exist 
for torus-like boundary conditions [57] and for a similar three-state QW j^ (but not in 
a truly honeycomb geometry) , the present is the first general treatment for the problem. 

The labeling of the sites in both formulations is indicated in Fig. 9 (b). Without 
loss of generality, we adopt the convention: in the y-direction the quantum number k is 
chosen such that if k is even, then the corresponding infinite row of sites have the bonds 
along north (up), southeast and southwest. On the other hand, for k odd, the bonds 
configuration at each site is south (down), northeast and northwest. 

An interesting aspect of the honeycomb lattice, consequence of how it imposes the 
relation between the quantum number a and directions, is that any sequence: 



(jo, ko)a -> (jl, ki)^ -^ (j2, k2)a ■■■ ^ (JN-I, 



u 



CAr_l 



a -^ iJN, k 



N 



(37) 
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(b) 



j-3j-2j-l j j+lj+2j+3 




k + 2: even 
k+1: odd 

k: even 
k-1: odd 



a = [j + 1]3 (7 = [j - 1]3 

Figure 9. (a) The honeycomb lattice and (b) the convention for the sites labehng. 
(c) For the coined formulation, the association between the quantum number a and 
directions. Here [x]^ denotes x mod 3, i.e., the remainder of x/3. 

i.e., a particular N steps evolution for which the value of a remains the same, is not 
naturally (i.e., necessarily) a ballistic-like trajectory. By natural we mean those cases 
(like regular lattices with an even coordination number, e.g., our square lattice), where 
we always can associate quantum states and directions such that a — )■ cr leads to 
evolution along straight lines. In such cases, the above dynamics, Eq. (!37|) . would 
be ballistic in the Hilbert lattice space regardless the specific a. As we are going to see 
for the construction adopted here, successive transitions cr — t- a yield a return to a same 
site (in a round trip) after visiting six sites. Hence, diffusion throughout the lattice 
implies transitions of the type o ^r o' ^ o . 



3.2.1. The coined model 

For the inner coin states |o"), we have three possibilities, namely, a = 0,1,2 (we 
also could call them 1,2,3, but the use of instead of 3 simplifies the notation). Due 
to the honeycomb particular topology, some care is necessary in defining the dynamics 
in terms of the coin variable. Indeed, states with a same a but at different sites does 
not always evolve to a same direction. So, to properly associate the a's with the system 
evolution under U, we consider the following prescription (which, however, is not the 
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k-2 



J - 3 j - 2 j - 1 j j + 1 j + 2 j + 3 



Figure 10. An illustration of the resulting association between coin quantum numbers 
and directions according to the rules in Fig. 9 (c). Here k is even and [j]^ = 0. 



only possiblqll). Let [x]n = x mod A^, i.e., [xJtv is the remainder of x/N . Then, for a 
given site (j, /c), the quantum numbers 0, 1, 2 are associated to directions in the lattice 
as indicated in Fig. 9 (c). An example of the resulting configuration is displayed in Fig. 
10 (in special, note from Fig. 10 that leaving from a certain site, say (j. A;), and always 
evolving to a same a, mandatorily will get back to (j, k) in exactly six steps). The rules 
in Fig. 9 (c) establish in an unique and self-consistent way - for the entire lattice - how 
single steps (j, k) — )■ (j', k') are determined by the values of a. 

The basis states are written as |(j, k)) ® |o"), with a column vector representation 
for \a) given by 

/o\ 



/l\ /o\ 

|0)= ,|1)= 1 

vv V) 

So, the coin operator at any site (j, k) is 



|2) 







(3^ 



/ ihk) {j,k) 



(jU,k) 







01 



c 



iJ,k) (j,k) 



10 



11 



V „{j,k) (j,k) 
\ '-2 '-21 



'-0 2 > 
'-12 
'-2 2 



(39) 



which we suppose to be an unitary matrix. 

Assuming the above construction, after a little lengthy but straightforward analyzes 
one finds that the one step time evolution operator for the coined formulation reads 

U.= (E^-®I^)(^|) X {j2\(j,k)){ij,k)\®C'^'''^), (40) 



vcr=0 



j,k 



where 



Sa\{j,k)) = \{f{j,k;a),g{j,k;a))), 

Sj\U,k))= \{f{j,k;M^)),9U,k;M^m, (41) 

I A full classification of equivalent constructions for the honeycomb lattice will appear elsewhere. 
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(j + l,k+l) 




(j,k-l) 



Figure 11. For the honeycomb lattice, the SQW basis states labeling rules. The a 
values here {aa, <Jb and CTc) are defined in Eq. (|44|) . The upper details summarize the 
convention in Fig. 9 (c) to label the CQW inner states. 



for (with sgn[x] = a;/|x|, if x 7^ 0, and sgn[0] = 0) 

/(j, k- a)=j + (-l)-b-+W^l3 sgn[a - [j + [k],],], 
gU, k-a)=k + {-If (1 - 2 sgn [\a - [j + [k]2];\] 
0fc(a) =[a-{-l)%. 

By construction 5*0- S^'^ = Sj^ S^ = I, leading to an unitary Uc- 
Finally, the corresponding projector operator is 

2 
^P^ = l(j,fc))((j,A;)|®5^|a)(a|. 

(7=0 



(42) 



(43) 



3.2.2. The scattering model 

For the scattering formulation, we consider the schematics in Fig 11. We denote all 
the basis state incoming to the site (j, fc) (continuous arrows in Fig. 11) by |cr, {j,k)). 
Since now the lattice coordination is 3, we need more than the usual ± to label the 
direction quantum number a. Therefore, we assume that a has the values 0,1,2, in 
analogy with the coined model. But in each bond we can ascribe only two of these 
possible values to the two states representing opposite directions. In this way, for the 
general situation represented in Fig. 11, we set the convention (which can be checked 
to be self-consistent along the whole lattice) 



C^a = [Jh, O-fe = [j + 1]3, CTc = [j - 1]3- 



(44) 



At each bond, for a state \a,{j,k)) incoming to {j,k) there is a corresponding 
outgoing state, which is itself an incoming state to another site (j', k') and denoted by 
|cr', (j', k')) (dashed arrows in Fig. 11). From the protocol used to label the cr's - Fig. 
11 and Eq. (jUj) - it is not difficult to show that a' = (pki^'), j' = f{j, k; 0^(0")), and k' = 
g{j, k; 0jt(cr)), for 0, / and g given in Eq. P2l) . For instance, in Fig. 11 for k even and 
[j]s = 0, ctq = and the other state in the same bond than |0, (j, k)) is |2, (j + 1, fc — 1)), 
as verified by direct inspection. Likewise, this follows directly from a' = 0fc-evcn(O) = 2, 
/ = /(j, k; 2)|fc-even,[i]3=o = j + 1 and k' = g{j, k; 2))|fc-even,[i]3=o = k - I. 
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Thus, the evolution operator is Ug = T + R, with 

T\(T,{j,k)) = ^t^^^'^^^^^\(f)k{a),{f{j,k;(f)k{a)),g{j,k;(f)k{a)))), 

R\a, (j, k)) = r^'g)^^ \M^), {f{j, k- </)fc(a)), (?(j, fc; </)fc(a)))), (45) 

where furthermore 

rt|a,(j,A;)) = Y. [tiT''^'^''^^^'^'^''^'^^^^^^^ 

i?t|a, (j, A;)) = [,(/0-M.W),.0-,fc;<^.W))]* |^^(^)^ (^(^. ^. 0,(^)), ^(^'^ fc; 0,(a)))). (46) 

Now, if we define 

/ r'^hk) p(i,fc) p(j,fc) \ 
/ ■•- 00 -*- 01 -'-02 \ 

rO-,fc) = rO>) po.fe) rO,fc) ^ (47) 

, rO'.fc) r(-?''^) r(-'''^) / 

\-'-20 -'-21 -'-22 / 

and identify (a, /3 = 0, 1, 2) 

^0.H,. = r^;'!)., t;,,, = r^^;') (for /? ^ 0,(a)), (48) 

one has that the unitarity of F in Eq. (H71) guarantees that t/^ is also unitary. 
Finally, the projector operator reads 

Pi^^^^ = \M^), {f{j, k- 0fc(a)), ^(j, k- 0fc(a))))((/(j, k- 0fc(a)), ^(j, k; 0fc(a))), 0fc(a)| 

+ |a,(j,A;))((j,A;),cr|. (49) 

3.2.3. Mapping the models 

First, we note that: (i) comparing the a labeling convention for the coined and 
scattering formulation in the schematics in Fig. 11 (see also Eq. (l44l) ): and (ii) taking 
into account how the states evolve according to the quantum number a in both models; 
it turns out that a direct one-to-one association between basis states is simply given by 

E\a,{j,k)) = \{j,k))®\a). (50) 

Second, if as done in Sec. 3.1, we set the coin matrix Eq. (139!) and the scattering 
coefficients in Eq. P7|) to be equal, then we find that the corresponding expressions for 
Us and f/c, with E given by Eq. (1501) . satisfies to Eq. (ITTj) . 

Lastly, to obtain the probabilities of one model by means of the other, we define 
the cross operators as in Eq. ( TTSj) . using the definitions in Eqs. ( H3l) . f H9l) and ( l50l) . 

3.2.4. Examples 

To illustrate the dynamics in a honeycomb topology, we again consider different 
coin (and equivalent scattering) matrices and calculate the QWs time evolutions. We 
analyze the following five different examples for the transition probabilities (assumed to 
be the same at all the lattice sites): 
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(a) 3 X 3 typical unbiased operator 

1 ( -' 

Cunba = -/= exp[-m/3] 

\ exp[— 7rz/3] exp 

(b) 3 X 3 biased (and real) operator 

/ 1 1-Vs 



exp[— 7r?/3] exp[— 7ri/3] 



-1 

-7ri/3] 



exp[— 7ri/3] 



(51) 



Q 



bia-j 



\ 






V3 



(c) 3 X 3 Discrete Hartley Transform (DHT) 

/2 

2 



C 



DHTs 



2^3 



V 



2 
1 + V3 
1- V3 






2 \ 

-1- v^ 
-1 + V3 / 



(52) 



(53) 



(d) 3 X 3 Grover 



C, 



G3 



1 

3 



-1 



-1 



2 
2 



(54) 



2 2 -1 y 

(e) 3 X 3 Discrete Fourier Transform (DFT) 



Cdft3 — —:= 



1 
%73 



exp[ 



-27rz/3] \ 



(55) 



exp[27ri/3] y 



/ exp[27rV3] 

1 
\ exp[— 27ri/3] 

Few comments about the matrices above are in order. The first case, Eq. flSTj) . is a 
typical (complex) unbiased operator since the coefficients have the same value for their 
modulus square, |co-'o-p = 1/3. Therefore, the resulting transition probabilities are equal 
to 1/3. On the other hand, the second, Eq. (152!) - here with all entries real numbers - 
is strongly biased because the transition probabilities are very different (but obviously 
summing up to 1). For instance, for the scattering formulation, fc-even and [j]^ = 
in Fig. 11 (cf. Eq. (iHD), we have |to,oP = 1/9 ^ 0.11, |ti,oP = 2(2 + V3)/9 ^ 0.83, 
^2,0^ = 2(2 — -\/3)/9 ~ 0.06. The discrete Hartley transform matrix [59j, Eq. ( |53|) . 
although not usually considered in QWs, is an interesting example to study due to its 
usefulness in signal processing [SU], furthermore always being real. The other two are 
just the 3x3 versions of the Grover and DFT. Finally, note there is not a 3 x 3 Hadamard 
matrix. 

Assuming as the initial states 

1^(0)), = |1, (0,0)), ^I^(O)), = |^(0))e = 1(0,0)) ® |1), (56) 

we evolve the models 20 time steps. The resulting probability patterns are displayed in 
Figs. 12-16. The graphs are organized exactly as done in the Sec. 3.1.4 (Figs. 5-8). 
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Figure 12. Quantum walks in tlie lioneycomb lattice considering the 3x3 typical 
unbiased transition coefficients, Eq. ([5T|) . The density plots represent the probabilities 
to be in the different states (as defined by the projector operators) after n = 20 
time steps for |vl'(0)) given in the main text. The results are for the: (a) coined, (b) 
scattering obtained from the coined, (c) scattering, and (d) coined obtained form the 
scattering, models. The graphs (a) and (d) and (b) and (c) are complete identical. 



Again, certain facts are clearly observed from Figs. 12-16. First, as it should 
be, we can obtain one evolution from the other by correct projections. Second, like 
in the square lattice, in the honeycomb the CQW probabilities are in general more 
sparse and somewhat smoother than those for the SQW. And third, it is interesting 
to notice the particular pattern in Fig. 13, with a tendency of three preferable 
directions of propagation along the lattice. This is a consequence of the probabilities 
bias resulting from the matrix Cbia-j, Eq- (l52l) . for which three transition elements 
Co2 = Cio = C21 = (1 + v^)/3 are considerable higher than the other six. 

Finally, comparing the Grover (DFT) probability amplitudes for the square and 
honeycomb topologies, respectively. Figs. 7 and 15 (Figs. 8 and 16) we can have an 
idea on the influence of a regular lattice coordination number to the QW dynamics. In 
fact, note that the |^E'(0)) used for the square lattice examples is equally "distributed" 
among the four a values basis states, Eq. f l36|) . which is not the case for the honeycomb 
lattice whose |\E'(0)) is written in terms of just one value of cr (from the three possible), 
Eq. (!56|) . Even then, we see that in the latter the resulting probabilities are more 
uniformly distributed, specially in the Grover case. This illustrates the intricate process 
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Figure 13. The same as in Fig. 12, but for the 3x3 biased real matrix, Eq. ([52 




Figure 14. The same as in Fig. 12, but for the 3 x 3 DHT, Eq. dSS]). 
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Figure 15. The same as in Fig. 12, but for the 3x3 Grover, Eq. ([M 




Figure 16. The same as in Fig. 12, but for the 3 x 3 DFT, Eq. §5 
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of quantum transitions at each lattice site as function of the number of bonds attached 
to it. 

4. Conclusion 

By means of detailed analysis of specific examples, in this work we have compared the 
different aspects associated to the construction of the discrete time coined and scattering 
quantum walks. Specifically, we have presented a complete formulation for the line, 
square and the involving honeycomb lattice topologies. In all the examples studied, we 
have clarified how to map one model version to the other. Also, we have illustrated 
the particularities in the probability pattern distributions along the graph structures, 
resulting from the different time evolution of each formulation. Finally, distinct coin 
(and corresponding scattering) matrices operators were considered. Besides the usual 
Hadamard, Grover and Discrete Fourier Transform, we have addressed the Discrete 
Hartley Transform as well as a few other cases. 

Our aim here was to show through concrete situations that in fact the CQWs 
and SQWs are closely related. Moreover, the use of a specific QW version in different 
applications may be more a matter of practicality (or even of taste, as recently stressed 
[H]), than due to prohibitive conceptual differences between the models. 

Regarding implementations, the issue may be a little subtle considering that usually 
distinct experimental architectures are proposed for CQWs and for SQWs (see the 
Introduction Sec). However, a given realization of a specific QW model conceivably 
could also be used to obtain results from the other model. But then, the experimental 
setup somehow should allow the effective construction of proper cross projectors (as 
defined in Sees. 2-3 and 2-4). The important point is that efforts in this direction 
certainly would pay off since an actual implementation of one version of QWs having 
proper built-in projections to the other version - thus "carrying" the features of both 
type of QWs - would be much more flexible and useful in applications. 
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